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Two-points problem for an evolutional first order equation in 

Banach space 
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Abstract 

Two-point nonlocal problem for the first order differential evolution equation with an operator co- 
efficient in a Banach space X is considered. An exponentially convergent algorithm is proposed and 
justified in assumption that the operator coefficient is strongly positive and some existence and unique- 
ness conditions are fulfilled. This algorithm leads to a system of linear equations that can be solved 
by fixed-point iteration. The algorithm provides exponentially convergence in time that in combination 
with fast algorithms on spatial variables can be efficient treating such problems. The efficiency of the 
proposed algorithms is demonstrated by numerical examples. 
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^ ! 1 Introduction 

The m-point initial (nonlocal) problem for a differential equation with the nonlocal condition u(io) + 
g{t\] . . . ; t p \ u) = uq and a given function g on a given point set P = {0 = to < t\ < ■ ■ ■ < t p } , is 
one of the important topics in the study of differential equations. Interest in such problems originates 
mainly from some physical problems with a control of the solution at P. For example, when the function 
■^j- ' g(ti] ■ ■ -]t p ]u) is linear, we will have the periodic problem u(io) = u{t\). Problems with nonlocal conditions 

arise in the theory of physics of plasma [19 , nuclear physics [TJ], mathematical chemistry [15], waveguides 
[T2] etc. Two-point problem is also useful for considering the finale value problem [22] . 

Differential equations with operator coefficients in some Hilbert or Banach space can be considered as 
meta-models for systems of partial or ordinary differential equations and are suitable for investigations using 
tools of the functional analysis (see e.g. [U US])- Nonlocal problems can also be considered within this 
framework [2J|3]. 

Discretization methods for differential equations in Banach and Hilbert spaces were intensively studied 
in the last decade (see e.g. El El EH ED HOI HD \M W\ and the references therein). Methods from 
[3 HI [9l HU [TTJ [211 [2S1 [27] possess an exponential convergence rate, i.e. the error estimate in an appropriate 
norm is of the type 0(e~ N ), a > with respect to a discretization parameter N — s> oo. For a given tolerance 
e such discretization provides optimal or nearly optimal computational complexity [7]. 

In the present paper we consider the problem 

*f + *(.«.)- AM, (L1) 

w(0) + au(l) = ip, 

where Ai{t) is a densely defined closed (unbounded) operator with the domain D(A{) independent of t in 
a Banach space X, ip is given vector and /i(i) is given vector- valued function, a£l We suppose that the 
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operator Ai (t) is strongly positive; i.e. there exists a positive constant Mr independent of t such that on 
the rays and outside a sector T,g = {z G C : < arg(z) < 6,6 G (0,7r/2)} the following resolvent estimate 
holds: 



(zI-Arit))- 1 ^ 



M R 

l + \z\ 



This assumption implies that there exists a positive constant c K such that ( see [B], p. 103) 

\\Af(t)e— Al ® || < c K s~ K , s > 0, k > 0. 
Our further assumption is that there exists a real positive uj such that 



-sA 1 {t)\ 



< e 



Vs, t G [0, 1] 



(1.2) 

(1.3) 
(1.4) 



(see [IS], Corollary 3.8, p. 12, for corresponding assumptions on Ai(t)). Let us also assume that the following 
conditions hold true 



||[Ai(t) - A 1 (s)}A-' f {t)\\ < L hl \t -s\ Vt, s, < 7 < 1, 



(1.5) 



||A5 , (i)Ap , (s)-J|| <L 7 |t-s| Vt, sG [0,1]. (1.6) 

We suppose also that 

/i(t)€C(0,l;I), (1.7) 

The aim of this paper is to construct an exponentially convergent approximation to the problem (jl.ll) 
for a differential equation with two-points nonlocal condition in abstract setting. The paper is organized as 
follows. In Section [2] we discuss the existence and uniqueness of the solution as well as its representation 
through input data. A numerical algorithm is presented in section [3] The main result of this section is 
theorem 13 .31 about the convergence rate of the proposed discretization. The next section|4]we represent some 
numerical example which confirm theoretical results from the previous sections. 



2 Existence and uniqueness of the solution 

It is well known, that for a = the problem (II. ip has a unique solution under the assumptions (|1.2p - (|1.7 
(se e.g. [HIES)- This solution can be write down as follows: 



pt pt 

i(t) = U(t,0)u{0)+ U(t,s)fi(s)ds = U(t,0)tp+ U(t,s)fi(s)ds 
Jo Jo 



(2.1) 



where U(t, s) is an evolution operator that corresponds to (11.11) for a = 0. 

Let us study conditions when there is unique solution for the two-points problem (jl.ll) . We have from 

423D 

u(l) = 17(l,0)u(0)+ / U(l,s)f 1 (s)ds. 

Jo 

Substituting this expression into the nonlocal condition we obtain 



u(0) = [1 + 0:17(1,0)] * 



tp-a l U(l,s)fi(s)ds 
o 



and for u(t) we have 



u(t) = U{t 1 0)[I + aU(l,0)Y 



tp — a 



U(l,s)f 1 (s)ds 



+ / U(t,s)f 1 (s)ds. 



(2.2) 



It is necessary to establish conditions on a for the existence of u(t). In fact, we have to explore when 
exists [/ + aU(l, 0)]~ . So, we obtain using estimate for U(t, s) (see e.g. [l8lfT3] h 

|[/ + al7(l,0)]- 1 || < [1 - \a\ FMU]- 1 < [1 - \a\M\- 1 < C, 

for small enough a (a < M~ x ). 



3 Numerical algorithm 

We use the approach developed in [10] and [25] to construct numerical method for solving problem (|1.1[) . 
First of all we change variable in (|1.1[) by t — >• -^ and for v(t) — u (^4^) we have 

f (— 1) + aw (1) = (p, 

where A(t) = \A X (±±i) , /(*) = i/j (±±*) , 

We choose a mesh w n = {£/,., k = 0, ...,n} of n + 1 various points on [— 1, 1] that are Chebyshev-Gauss- 
Lobatto nodes tk = cos (^^tt) , and set Tk = tk — th-x- Let 

3(t)=Afc = A(t fc ),tG(tfc_i,t fc ], fe = M, 

A = A(-1). (3 ' 2) 

Let us rewrite the problem (|3.1[) in the equivalent form 

^+3(t)« = p(t)-^(t)Mt) + /(*), te(-i,i) 

u(— 1) = ip — av{l). 

Note, that now all operators on the left side of these equations are constant on each subinterval and 
piece- wise constant on the whole interval [—1,1]. 

On each subinterval we can write down the equivalent to (|3.3[) integral equation 

v(t) =e- Ak( - t - tk - l) v(t k - 1 ) + f e- Ak(t - s ^[A k ~A(t)]v(s)ds+ f e - A ^- s ) f(s)ds, 

Jtk-i J U-x (3-4) 

t e [tk-i,tk], k = 2,n, 

v(t) = e- A ^ t+1 ^ [ip - av(l)} + f e- Alit - 8) [A l -A(t)]v(s)ds+ [ e~ A ^- s) f{s)ds, i € [-Mi]. (3.5) 



Let 

n 

P n (t;v) = P n v = Y,v(jk j )L j , n (t), (3.6) 

3=0 

be the interpolation polynomial for v(t) on the mesh u) n , x = (xq, ...,x n ),Xi S X given vector and 

P n (t;y) = P n x = ^XjLj, n (t) (3.7) 

j=o 

the polynomial that interpolates x, where 

£,•»(*) = "3 y " (s)(1 ~ 52) , j = 0,...,n 

are the Lagrange fundamental polynomials. Substituting P n (s;x) for v(s), Xk for v(tk) and then setting 
t = tk in (|3.4[) we arrive at the following system of linear equations with respect to the unknown Xk ■ 

x + ax n = (p, 

— Auti v~^ i i i — (3-8) 

x k = e k x k -i + ^akjXj +4>k, k = l,n, 

3=0 



which represents our algorithm. Here we use the notations 



a>kj = / e 
r-tfe 



Ok — I c 



-^('*-*)[A fc -A( s )]% n ( s )d S) 



(3.9) 



and suppose that we have an algorithm to compute these coefficients. 

For the error z = (zx, ..., z n ), with Zk = v(tk) — Xk we have the relations 

z + az„ = 0, 

n 

z fe = e~ AfcTA: z fc _i + )]akjZj +4>k, k = l,n, 



j'=o 



where 



^fc = / 



■ A *(**-')[Afc-A(*)][t;(*)-P„(»;t;)]d8, fe = l,n, 



(3.10) 



(3.11) 



In order to represent algorithm (|3.8[) in a block-matrix form we introduce the matrix 



S = 



( I 

-ax I 

-o-2 J 

\ 



cmto\ 





-O n I ) 



(3.12) 



where a = A^A„ 7 , cr*,. = e j4fcrfc AjM fe j' :L , fc = l,n, the matrix 5 = {a fcj }£ J=0 with a k j = A^afejA^ 
fe = 1, n, j = 0, n, and <Soj = 0, j = 0, n, the vectors 



AJxi 



( A^\ 
Alfa 



(Alzo\ 

AJ Zl 



i/> 



\Alx n j \AIK ) \Alz n j 

It is easy to check that for the (left) inverse 

S- 1 = S (Rx - R 2 ) , 



( o \ 

AJfa 



\*M>n J 



(3.13) 



(3.14) 



where 



Ri = 



S = (I + aa Q ai ...a n ) , 

( I o 

ox I 

<72<Tl 2 



\a n ■ ■ ■ ax o n - ■ -a 2 

/0 a n . . . o 2 o n . . .<j 3 
ax(Tn • ■ • CT 3 



R 2 = as 





\0 



0\ 



O n ij 

O n 
OxO n 



I \ 

Ox 



<7 n _i...<7l 

/ 



Remark 3.1 Using results of JllL ^llffl one can get a parallel and sparse approximations with an exponential 
convergence rate of the operator exponentials contained in S^ 1 and as a consequence a parallel and sparse 
approximation of S^ 1 . 



We multiply the equations in (|3.8[) and the equation in p.lOj) by AJ, k — 0, n and obtain 
A^xq + aAlx n = A^ip, 

n 

Alxk = e- AhTk Ajx k -i + Y^ & kjA] Xj + Alfa, k = T~R, 

3=0 



A^z + aA2z n = 0, 

n 

A\z k = e' AkTk Apk-i + J2 a kjA]zj + Altp k , k = I~n, 



(3.15) 



(3.16) 



Then the systems (|3.15|) . (|3.16[) can be written down in the matrix form using notations (|3.12|l . (|3. 13|) as 

Sx = Bx + <b, , 

(3.17) 

Next, for a vector v — («i, i>2, •••, v n ) T and a block operator matrix A — {a,j}™ =1 we introduce a vector 
norm 



and the consistent matrix norm 



\\\ v \\\ = \\\ v \\\i = max \\v k \\, 

l<k<n 

it 

IIHII = IIHI|i= m^y^llaijll- 

Ki<n *• — ' 
~ " 3=1 

Due to (ESI) we have IH^A^JI = p^"^ - / + /||| < 1 + L 7 r fc , ||<t || = H4^n 7 ll < 1 + L^T. In our 
case T = 2. So, we have the following, using these estimates 

|K || = Ue-^MZ^J < e-^||^Vill < ^ Tk 0- + V*) , 

\\S\\ = || (1+ aa Q a x . ..a^ 1 || < (1 - \a\ ||<r || |ki|| ||<7 2 || ■ • ■ KID" 1 

< (1 - \a\ (1 + 2L 7 ) e-" ri (1 + L 7 n) e^ 2 (1 + L 7 r 2 ) . . . e -" T " (1 + L 7 r n )) _1 

< ( 1 - |a| (1 + 2L 7 ) e- 2w fl + ^-) ) < (l - |a| (1 + 2L 7 ) e - 2w e 2 ^) _1 < c, 



for a small enough. 

In order to estimate the norm of matrix S we have to estimate the norms of matrices R\ , R2 ■ In 110] it 
was proved that for matrix similar to R\ the estimate |||i?i||| < en holds true. Let us estimate the norm of 
matrix i? 2 . 

|||i?2||| < (1 + 2c) (1 + e- WT (l + ct) + • ■ • + [e- WT (l + ct)] 11 - 1 ) 

/ fix rr V" — 1 \ p 2c 
< (1 + 2c) 1 + (1 + ct) + • • • + (1 + ct)' 1 - 1 < y ' < (1 + 2c) < en. 

\ CT J CT 

Using these estimates we obtain that 

IHS^III <cn. (3.18) 

It was proved an estimate for the matrix B in [10] : 

\\\B\\\<cn f - 2 \n(n). (3.19) 

So we can formulate the following assertion 
Lemma 3.2 Let the assumptions (|1.2p - (|1.6[) are fulfilled then the estimates (|3.18|) . (|3.19p hold true. 



Using (|3.17p we have 

L , (3-20) 

where _B is a diagonal matrix with unit operators / on diagonal. Using lemma 13.21 we obtain that 

IIIS^BIH < en 7-1 ln(n) -» 0, n^oo. (3.21) 

It means that for n large enough there exists the matrix \E — S~ 1 B] and 



[e-s^bY 1 



<c. 



Consequently we obtain the following stability estimates from p.20p using lemma [32 

p|||<cn|||0|||, 



-:i I : n,y ■" ! ' > " 2i; 



Let n„ be the set of all polynomials in t with vector coefficients of degree less or equal then n. In complete 
analogy with [U [531 121] the following Lebesgue inequality for vector- valued functions can be proved 

\\u(t)-P n (t;u)\\ C [-i,i]= max \\u(t) - P n {t;u)\\ < (1 + A„)E n (u), (3.23) 

te[-i,i] 

with the error of the best approximation of u by polynomials of degree not greater then n 

EJu)= inf max \\u(t) - p(t)\\. (3.24) 

pen„ {£[-1,1] 

Now, we can go over to the main result of this section. 

Theorem 3.3 Let the assumptions of Lemma \3.2\ with 7 < 1 hold, then there exists a positive constant c 
such that 

1. For n large enough it holds 

PHI Ken 1 - 1 - Inn • E n (A&v), (3.25) 

where v is the solution of (|3.1[) ; 

2. The first equation in (|3.17[) can be written in the form 

x = S^Bx + S^cj), (3.26) 

and can be solved by the fixed point iteration 

-(fe+i) = s -i B %(k) + s -i^ fe = 0, 1, ...; £(°) - arbitrary, (3.27) 

with the convergence rate of an geometrical progression with the denominator q < en 7 ™ 1 ln(n) < 1 for 
n large enough. 

Proof. For z we have the second estimate in (|3.22[) . The norm of the first summand on the right hand side 
of this inequality can be estimated in the following way 



max 

Kk<n 



ds 



{Alc- A ^-^[A k -A(s)}AP(AlA-')(A2v( s )-P n (s;A2v))} 

1 

(' '' k 

< c max / \t k -,s|" 7 |t fc ~s\ \\A^v(s) - P n {s; A^v)\\ds 

< cr^l \\Alu{s) - P„(-;4Ju)|| c[ _ M] < cr^- 7 (l + k n )E n {A»- 
So, we obtain 

III^HI < en 1 ' 2 -\nn-E n {Alu), (3.28) 

Now, the first assertion of the theorem follows from (J3.22I) , p.28[) . The second one follows from f|3.1T[) 
and P~2T1) . ■ 



4 Examples 

Let us consider the following problem 

du(x,t) d 2 u(x,t) 



dt dx 2 

u(0,t) =u{l,t) = 0, 
u(x, —1) + au(x, 1) = <p{x), 



+ q(x,t)u(x,t) = f(x,t), 



(4.1) 



with f(x,t) = e _7r2(1+t) sin(7rx)(l + i), a = 0.5, ^(») = [l + 0.5e~ 27r2 ") sin(7ra;), q(x,t) = 1+t. Then, the 

solution of this problem is u(x, t) = e _7r ( 1+t ) sin(7ra;). 

The problem (14.111 can be write down in the form (|3.1[) where the operator A(t) is defined by 



D(A(t)) = D(A) ={ue iJ 2 (0, 1) : v(0) = 0, w(l) = 0}, 



(4.2) 



Coefficients of the system (|3.15|l were calculated by using the Fourier series expansion. The results of 
calculation are presented in tables confirm our theory above. 



Point t 


£ 


-1 


0.00005276 


-0.70710678 


0.00097645 





0.00063440 


0.70710678 


0.00029592 


1 


0.00010552 



Table 4.1: The error in the case n = 4, x = 0.5 



Point t 


e 


-1 


8.12568908Ee-7 


-0.86602540 


0.00010146 


-0.5 


0.00030932 





0.00022136 


0.5 


0.00013419 


0.86602540 


0.00007182 


1 


0.00000162 



Table 4.2: The error in the case n = 6, x = 0.5 
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